Fractional quantum Hall states at zero magnetic field 
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We present a simple prescription to flatten isolated Bloch bands with a nonzero Chern number. 
We first show that approximate flattening of bands with a nonzero Chern number is possible by 
tuning ratios of nearest-neighbor and next-nearest-neighbor hoppings in the Haldane model and, 
similarly, in the chiral-vr-flux square lattice model. Then we show that perfect flattening can be 
attained with further range hoppings that decrease exponentially with distance. Finally, we add 
interactions to the model and present exact diagonalization results for a small system at 1/3 filling 
that support (i) the existence of a spectral gap, (ii) that the ground state is a topological state, and 
(iii) that the Hall conductance is quantized. 



In a seminal paper, Haldane [1 has shown that nonin- 
teracting electrons hopping on a honeycomb lattice can 
exhibit the integer quantum Hall effect (IQHE) with- 
out the Landau levels induced by a uniform magnetic 
field, provided the system breaks time-reversal symme- 
try (TRS). In that model, electrons hop with a real- 
valued uniform nearest-neighbor (NN) amplitude of mag- 
nitude that preserves TRS, as well as complex- valued 
next-nearest-neighbor (NNN) amplitudes with the uni- 
form magnitude ^2 that break TRS. A non- vanishing ^2 
generically opens a band gap at the Fermi-Dirac points 
of graphene (half- filling). This band gap results in the 
Chern numbers taking opposite values of magnitude 1 
on the upper and lower bands. Consequently, the model 
exhibits an IQHE at half-filling. 

Given the fact that a band insulator can support the 
IQHE without a magnetic field, a natural question that 
we address in this Letter is whether a fractional quantum 
Hall effect (FQHE) is also possible in an interacting lat- 
tice model without a magnetic field. For the usual FQHE 
in an uniform magnetic field, all Landau levels share the 
same Chern number, ±1 depending on the orientation of 
the uniform magnetic field. Moreover, in the absence of 
disorder, all Landau levels are flat (i.e., dispersionless) 
and thus can accommodate, when partially filled, an ex- 
ponentially large number of Slater determinants, from 
which incompressible liquids are selected by interactions 
at some special filling fractions. Haldane's model fulfills 
the first ingredient for the FQHE: nonvanishing Chern 
numbers for the single-particle Bloch bands. We are go- 
ing to construct two-dimensional lattice models without 
magnetic fields that also satisfy the second ingredient for 
the FQHE: band flattening. 

There is a long history of flatband models. They have 
been studied since the 1970s in amorphous semiconduc- 
tors [2H3], and understood using projection operators [5]. 
More recently they have been studied on kagome, honey- 
comb, and square lattices j6HTT]. In Ref. [lOj flatbands 
were isolated by gaps, and the question of whether it is 
possible to have a flatband with nonzero Chern number 



was raised. We shall answer this question affirmatively. 
We then add interactions and show evidence that the 
many-body state is a topological state with fractional 
Hall conductance at 1/3 filling. 

Our starting point is two-dimensional local lattice 
models describing the hopping of spinless fermions. In 
the spirit of Haldane's model, we restrict the lattice mod- 
els to those with only two Bloch bands and enforce lo- 
cality by only allowing NN and NNN hoppings. We will 
show that, by varying the ratio of the NNN to NN hop- 
pings, we can deform the bands to make them flatter. 
The characteristic measure for the flatness of a Bloch 
band is here the ratio of the bandwith to the band gap. 
We then show that this criterion for flatness can be sat- 
urated to the ideal limit of zero for the valence band by 
including arbitrary range hoppings. However, the flat- 
tened single-particle Hamiltonian still preserves locality 
in the sense that the hopping amplitudes decrease expo- 
nentially with the distance between any two lattice sites. 

Consider the noninteracting two-band Bloch Hamilto- 
nian of the generic form 

^0 •= ^i^k^k^ ^k '= ^cfcCTo ^Bj^ a. (la) 

keBZ 

Here, BZ stands for the Brillouin zone, = 
(^Ia'^Ib)' where ^ creates a Bloch states on sub- 
lattice s = A, B, and the 2x2 matrices ctq and a are the 
identity matrix and the three Pauli matrices acting on 
the sublattice indices. If we define 

'= rs^^ tan^fc := cos l9fc := Boj^, (lb) 

we can write the eigenvalues of Hamiltonian TL^, ass_^_ = 
Bq \Bi^ I and for the corresponding orthonormal eigen- 
vectors 

_ / e-''^'e/2 cos %\ _ ( e-'"^^!'^ sin % \ 
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Two examples of Hamiltonians of the form (la) are the 
fohowing. 

Example 1: The honeycomb lattice. — We introduce 
the vectors a\ = (0,-1), a\ = (^3/2, 1/2), a\ = 
(— >/3/2, 1/2) connecting NN and the vectors b\ = 



a 



3' 



^2 



a*^ — b*^ = a\ — 02 connecting NNN from the 
honeycomb lattice depicted in Fig.jlja). We denote with 
k a wave vector from the BZ of the reciprocal lattice dual 
to the triangular lattice spanned by and 63, say. The 
model is then defined by the Bloch Hamiltonian [T] 



(2a) 
(2b) 




ti COS k • 
ti sin k • 



=1 \ — 2^2 sin<l>sin/c • 6^ 



where > and ^2 > are NN and NNN hop- 
ping amplitudes, respectively, and the real numbers ±$ 
are the magnetic fiuxes penetrating the two halves of 
the hexagonal unit cell. For t^ ^ ^2, the gap A = 
min^e^ ^ — max^6:_ ^ is proportional to ^2- The width 
of the lower band is 5_ = max^£:_ ^ — min^£_ ^. The 
fiatness ratio ^_ / A is extremal for the choice cos <l> = 
^1/(4^2) = 3^^3/43, yielding an almost fiat lower band 
with (5_/A = 1/7 [see Fig.[^c)]. 

Example 2: The square lattice. — We introduce the 
vectors = (l/V^, 1/V^) and = (-1/^2, 1/V^) 
connecting NNN from the square lattice as depicted in 
Fig. [l|b). We denote with k^ = (k^^ky) a wave vector 
from the BZ of the reciprocal lattice dual to the square 
lattice spanned by x and y. The model is then defined 
by the Bloch Hamiltonian [12] 



-i7r/4 



1 + e+' 



^3,fc := 2*2 (cosfc^ -cosfcy 



(3a) 
(3b) 
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where t-^ > and ^2 > are NN and NNN hopping ampli- 
tudes, respectively. The fiatness ratio S_/A is extremal 
for the choice ^1/^2 = a/2, yielding two almost fiat bands 
with S_/A ^ 1/5 [see Fig.jljd)]. 

The Chern numbers for the bands labeled by ± in 
Eq. (lie]) are given by 
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d, cos 8{k) d,m . (4) 
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They have opposite signs if nonzero. All the information 
about the topology of the Bloch bands of a gaped sys- 
tem is encoded in the single-particle wave functions. For 
example, the Chern numbers depend solely on the eigen- 
functions. Haldane's model ([2| and the chiral-7r-fiux ([3| 
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FIG. 1. (Color online) (a) Unit cell of Haldane's model on 
the honeycomb lattice: the NN hopping amplitudes ti are real 
(solid lines) and the NNN hopping amplitudes are t2e^^^^^^° 
in the direction of the arrow (dotted lines). The flux 3$ 
and — $ penetrate the dark shaded region and each of the 
light shaded regions, respectively. For $ = 7r/3, the model 
is gauge equivalent to having one flux quantum per unit cell, 
(b) The chiral-TT-flux on the square lattice, where the unit cell 
corresponds to the shaded area. The NN hopping amplitudes 
are tie^^/^ in the direction of the arrow (solid lines) and the 
NNN hopping amplitudes are t2 and — 12 along the dashed and 
dotted lines, respectively, (c) The band structure of Haldane's 
model for cos$ = ti/(4t2) = 3^3/43 with the flatness ratio 
1/7. (d) The band structure of the chiral-vr-flux for ti/t2 = 
with the flatness ratio 1/5. The lower bands can be made 
exactly flat by adding longer range hoppings. 



are topologically equivalent in the sense that both have 
two bands with Chern numbers ±1. 

To enhance the effect of interactions, highly degener- 
ate (i.e., fiat) bands are desirable. It is always possible 
to deform the Bloch Hamiltonian (la) so as to have one 
fiatband with the energy —1, say, while preserving the 



eigenspinors x± k (P^- Indeed, this is achieved by turn- 
ing the Bloch HamrHonian (la) into 



flat 



(5) 
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Note that whenever 5q ^ =0, the Hamiltonian (la) has 
the spectral symmetry £^ ^ = —£_ ^ so that botn bands 

of = Bj^ • a are completely fiat. This spectral sym- 
metry applies to the chiral-7r-fiux ([3| but not to Haldane's 
model ([2| unless ^ = ±7r/2. 

Generically, follows from a lattice model for which 
the hopping amplitudes are nonvanishing for arbitrary 
large separations. If, however, the hopping amplitudes 
decrease sufficiently fast with the separation, locality is 
preserved. To estimate the decay of the hopping ampli- 
tudes with the separation between any two sites in the 
ffattened chiral-7r-ffux model ([3|, we calculate the decay 



3 



of the Fourier coefficients ^, of 



Here, we introduce 



— 1= ^ A^,n/Cosn/c+cosn'/c_, (6) J ^^'^k ^-,k[^k ^ [xl^k^kX-,k) 



where = k^^ky. Because 

^±,fc = (2^? - ^2) (2 + cos k^ + cos k_) 
+ ^2 (3 + COS kj^ COS k_) , 

it is sufficient to consider the Fourier coefficients A^ of 



(9b) 



1 ^ ^ 

= A^cosn/c, for — 1 



Vl + a cos k 



< a < 1. (8) 



n=0 



In the hmit n ^ \^ one finds that A^ decays exponen- 
tiahy with n. We conclude that for any fixed n, the coef- 
ficients A^ ^/ decay exponentially with r\! and vice versa 
for n, iff \a\ < 1. Flattening the energy bands preserves 
the locality of the chiral-7r-fiux model ([3|. 

The fact that we have engineered single-particle wave- 
functions in a fiat Bloch band that support a Chern num- 
ber ±1 is one step in mimicking the FQHE. However, be- 
cause of lattice effects, it is not a given that interactions 
lead to a many-body ground state supporting the FQHE. 
In fact, even when there is a uniform magnetic field, the 
combination of lattice effects and interactions is not well 
understood upon increasing the magnetic fiux threading 
the elementary lattice unit cell. In Refs. [13", TT, for ex- 
ample, the possibility of a FQHE induced by interactions 
for the Hofstadter problem, NN hopping with a uniform 
fiux threading each elementary plaquette of the square 
lattice, was studied numerically. While the overlap be- 
tween the Laughlin states on the torus and the lattice 
many-body ground states was close to unity when the 
plaquette fiux is much smaller than the fiux quantum, 
this overlap decreases when the plaquette fiux becomes 
of the order of one quarter of the fiux quantum. It is 
thus imperative to study how interactions lift the macro- 
scopic degeneracy of a fractionally filled fiat Bloch band 
and whether a gapped topological ground state emerges. 

Two distinctive properties of such a ground state |^) 
at filling fraction u (where is an odd integer) and 
with periodic boundary conditions (toroidal geometry) 
are (i) the i^~^-fold topological degeneracy of the ground 
state manifold and (ii) the quantization of the Hall 
conductance a^y in units of ve^ /h. The Hall conduc- 
related to 



tance is related to the Chern-number C 
body ground state |^) as a^y = Ce'^/h 
tionally, the Chern-number is evaluated using 



of the many- 
15 . Conven- 
twisted 



boundary conditions (r + N^x]"^^) = e^^^(r|^^) and 
{r + NyVl^^) = e'^y{r\^^), where 7^ = (7^,7^;) is the 
twisting angle and x Ny the number of unit cells. The 
Chern number is then given by [16] 



C 



1 
2^1 



J dSV^A(^^|V^|^^). 



(9a) 



as a second way to calculate the Chern number, where 
n_ = (^|c^ _Cj^ -1^) is the occupation number of the 
single-particle Bloch state in the lower (— ) band with 
wave vector k evaluated in the many-body ground state. 
We can show that both formulas are equivalent, i.e., 
C = C. To this end, one expands the many-body wave 
function |^) in a sum over Slater determinants and ap- 
plies a gauge transformation to the single-particle states 
to remove the twist in the boundary conditions. 

We close this Letter with an exact diagonalization 
study to show the existence of a gapped topological 
ground state for the chiral 7r-fiux phase ([3| in the pres- 
ence of interactions. We consider an interaction defined 
by the repulsive two-body NN potential V^j according to 



(ij) 



V>0. (10) 



Te[0,27r] = 



Directed NN bonds of the square lattice A = A U B made 
of the open and filled circles of Fig.jljb) are here denoted 
by (u)? while is the occupation number on the site 
i e A. 

We also drive the model trough a topological phase 
transition to establish that a gapped topological many- 
body ground state results from the topological nature of 
the model. To this end, we add a sublattice-staggered 
chemical potential 4/ig to the single-particle Hamilto- 
nian ([T]) by replacing B^ j^ -\- 4/ig in Eq. ([3|. 

Then, the two noninteracting bands have a Chern num- 
ber ±1 for 1^2/^8! ^ 1 ^ vanishing Chern number 
for \t2/Ps\ < 1- The topological phase transition at 
1^2/Msl = 1 forces the single-particle spectral gap to close 
and the fiattening of the bands is ill-defined at that point 
[a = 1 in Eq. 

For a 3 X 6 sublattice A, we find a unique ground state 
that is separated by a gap of the order of the interac- 
tion strength V at filling u = 1/3 [see Fig. [2|a)]. This 
state loses its clear separation in energy from the other 
states at the topological phase transition \t2/Ps\ — ^ 
another gapped ground state is obtained for 1^2 //igl < 1. 
Based on (i) the topological degeneracy and (ii) the quan- 
tized Hall conductance, we will now argue that the first 
state is a topological many-body state while the latter 
state is topological trivial. 

(i) Because of translational invariance, the Hamilto- 
nian does not couple states with different center of mass 
momenta Q := /c^ + . . . + /c^, where k^^ i = 1, • • • , 
are the single-particle momenta of an A^-particle state. 
At 1/3 filling of the 3x6 sublattice A, the particle num- 
ber A/" = 6 is commensurate with the lattice dimensions 
and all three topological states have the same Q. As a 
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FIG. 2. (Color online) (a) The lowest eigenvalues of H^""^ + 
i^int foi" the chiral vr-flux phase obtained from exact diago- 
nalization for 6 particles on a 3 x 6 sublattice A (1/3 fill- 
ing), normalized by the bandwidth E^. The parameters t2 
and /is of Hq^^ interpolate between topological (|t2| > |Ms|) 
and nontopological (|t2| < |/is|) single-particle bands. Here, 
g (2/7r) arctan |t2//is| and the energies are measured rela- 
tive to the energy of the single-particle band, (b) The lowest 
eigenvalues in the center of mass momentum sector of the 
ground state of Hq^^ + H^^^ with twisted boundary conditions 
as a function of the twisting angle 7^, for ji^ = Q, t2 = ti/v^. 
The level crossings indicate the topological nontrivial nature 
of the three lowest states, (c) Same as (b), but for ji^ — ti/v^, 
t2 = 0. The ground state is topologically trivial. 



consequence, their topological degeneracy is lifted and a 
unique ground state appears. We can now use twisted 
boundary conditions to probe the topological nature of 
the ground state: varying 7^ between and 27r is equiv- 
alent to the adiabatic insertion of a flux quantum in the 
system. During this process, a topological ground state 
with (7 = 1/3 should undergo two level crossings with the 
other two gapped topological states [17 . Indeed, we find 
these level crossings for the gapped ground state when the 
model has topological single-particle bands [Fig. |2|b)], 
whereas no level crossings are found otherwise [Fig.|2[c)]. 

(ii) We have also calculated the Chern number of the 
gapped ground state for /ig = 0, t 

and 



2 — ti/\/2 with the 
l9bl). We find C = 



two equivalent formulas (9a 
0.29 and C = 0.30 and attribute the deviations from 
C = 1/3 to the limitations of the small system size. For 
the topological trivial model with /ig = ti/V2^ ^2 = 0, 
we find that C and C vanish to a precision of 10~^ and 
10~^, respectively. 

In summary, we have proposed a simple recipe to de- 



form any noninteracting lattice model so as to obtain 
fiatbands, while preserving locality. We flattened the 
bands of the chiral 7r-flux phase and then lifted the result- 
ing macroscopic ground state degeneracy with repulsive 
interactions. Via exact diagonalization, we have shown 
that a FQH-like topological ground state is obtained at 
1/3 filling. This ground state, that is not well described 
by Laughlin-type wavefunctions, will be further studied 
in future works. 
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Note added. — Recently, we became aware of Refs. [I8l 
IT9) in which similar topological flatband models are dis- 
cussed. Subsequently, Ref. [20] appeared with exact diag- 
onalization results that are consistent with our findings. 
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